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We present time- and angle-resolved photoemission spectroscopy measurements on the charge 
density wave system CeTea. Optical excitation transiently populates the unoccupied band structure 
and reveals a gap size of 2 A = 0.59 eV. The occupied Te-5p band dispersion is coherently modified 
by three modes at = 2.2 THz, Q 2 = 2.7 THz and Q 3 = 3 THz. All three modes lead to small 
rigid energy shifts whereas A is only affected by and Their spatial polarization is analyzed 
by fits of a transient model dispersion and DFT frozen phonon calculations. We conclude that the 
modes and D 2 result from in-plane ionic lattice motions, which modulate the charge order, and 
that D 3 originates from a generic out-of-plane Aig phonon. We thereby demonstrate how the rich 
information from trARPES allows identification of collective modes and their spatial polarization, 
which explains the mode-dependent coupling to charge order. 

PACS numbers: 78.47.J-, 71.45.Lr, 73.20.Mf, 63.20.kd 


Interaction of multiple degrees of freedom in quan¬ 
tum materials can result in ordering phenomena that 
give rise to complex phase diagrams with competing or 
intertwined orders [l|. Such broken-symmetry ground 
states host characteristic low-energy excitations, which 
are intimately related to the interactions responsible for 
ordering. Ultrafast non-equilibrium methods are a new 
and active field for the investigation of broken-symmetry 
phases as they grant direct access to collective modes in 
the time-domain [113 and can provide strong evidence 
for the ordering mechanism operating in a particular ma¬ 
terial [1,0. 

Gharge density wave (GDW) formation is one of the 
most basic and widely studied cases for long-range order¬ 
ing when only electronic and lattice degrees of freedom 
are interacting 0. Below a critical temperature Tcdw 
a spontaneous lattice distortion Az appears, generally 
originating via the complex interplay of the electronic 
structure and lattice. This results in a spatial charge 
modulation at a wave vector q (spatial period 27r/g) and 
opening of a bandgap 2A in the band structure, which 
minimizes the system’s total energy by reducing the ki¬ 
netic energy of the electronic system. One of the charac¬ 
teristic collective excitations in such Peierls-like GDWs 
is a gapped amplitude mode (AM), which modulates the 
magnitude A(t) of the order parameter 0. 


Quasi-2D rare-earth tritellurides RTes (i7=Y, La-Sm, 
Gd-Tm) form an incommensurate GDW along the crys¬ 
tallographic c-axis and have been studied extensively [111 - 
0 . At low temperatures, the heaviest members of the 
RTes series (R=Dy-Tm) additionally order with a second 
GDW perpendicular to the first GDW 000 . Angle- 


resolved photoemission spectroscopy (ARPES) reveals 
the Fermi surface (FS) top olog y and the momentum- 
dependent gap structure |2 Ch24|. 


Gollective modes in RTes were studied with Raman 
spectrosco Dv [iqI . 25, 20, time-resolved optical reflectiv- 
jiy ^ and time-resolved ARPES (trARPES) 

0, [ 29 I, [30|. Yet even for this thoroughly investigated 
GDW system unambiguous identification of an AM is 
non-trivial. Lowering the translational symmetry in the 
GDW phase leads to the renormalization of the normal 
state optical and acoustic phonon branches at g [25- 
27|. Ubiquitously observed generic phonon modes add 
to the difficulty of mode identification. Direct com¬ 
parison of modes observed in the GDW phase to cal¬ 
culated phonon spectra is hampered by the incommen¬ 
surable GDW in RTes. These challenges are present 
in GDW materials in general and have led to conflict¬ 
ing reports regarding the AM in RTes in particular. In 
DyTes, for instance, a 1.75 THz mode is observed in both 
temperature-dependent transient reflectivity [13] , Raman 
spectroscopy [20 and a trARPES experiment using a 
pump-pump-probe scheme [ 10 . Optical measurements 
attribute this mode to a phonon in the lower symmetry 
GDW state that is coupled to the AM, whereas trARPES 
[10 assigns it to a second AM. 


In this paper, we report a trARPES study of the col¬ 
lective response of GeTes in the GDW phase. The unoc¬ 
cupied band structure is transiently populated by optical 
transitions, which reveals a gap size of 2A = 0.59(2) eV. 
The occupied Te-5p band dispersion near the GDW gap is 
coherently modified by three modes and analyzed within 
a transient tight-binding (TB) model. Momentum de- 
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pendent analysis combined with DFT frozen phonon cal¬ 
culations reveals different coupling to the band structure 
for in-plane AMs and out-of-plane optical phonons. 

Small rigid shifts of the band dispersion Eq are in¬ 
troduced by all three modes, whereas much larger pe¬ 
riodic modulations of A are driven by lattice motions 
fli = 2.2 THz and ^2 = 2.7 THz. We explain this ob¬ 
servation by concluding that both modes and Q 2 re¬ 
sult from in-plane ionic motions which are coupled to 
the CDW and hence drive oscillations of A as AMs. The 
Qs = 3 THz mode originates from an out-of-plane optical 
Aig phonon that does not couple to the CDW. Our work 
demonstrates how careful evaluation of the rich informa¬ 
tion from trARPES allows unambiguous identification of 
AMs in charge-ordered materials. 



FIG. 1. (a) Structure of the Te-planes where the CDW forms 

along the c-axis with lattice displacement Az indicated. The 
FS is described in a TB model, with the overlap integrals tpa- 
along and tpn perpendicular to Te-5p orbitals, (b) Cartoon of 
the trARPES experiment on CeTes in the CDW region of the 
FS {dotted lines). The blue plane illustrates an energy versus 
momentum cut with Te-5p bands and a pump pulse exciting 
electrons into unoccupied states. The blue line indicates the 
momentum cut in the experiment across the gapped region 
{dotted lines) of the calculated 2D FS, green arrow denotes 
the CDW wave vector q. trARPES spectra for selected pump- 
probe delays, 330 fs before the arrival of the pump pulse (c) 
and during the presence of the pump pulse (d). Intensities 
are rescaled exponentially as function of energy for enhanced 
visibility of the transiently occupied features above Ep. (e) 
TB fit {blue markers) according to Eqs. ([T|) and © to the 
extracted band dispersion {blaek markers) at —25 fs yields 
2 A = 0.59(2) eV. Size of the blue markers is proportional to 
the calculated spectral weight of the TB bands [^. Brown 
lines indicate the calculated metallic bare band dispersion. 


Our trARPES setup 
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consists of an amplified 


laser system operating at 310 kHz repetition rate that 
drives an optical parametric amplifier to optically ex¬ 
cite the sample with ^ 50 fs pump pulses at a pho¬ 
ton energy of 1.03 eV. Electrons are photoemitted with 
^ 150 fs probe pulses at 6 eV as function of pump-probe 
delay and collected with 50 fs delay steps in a hemi¬ 
spherical electron analyzer, yielding a total energy res¬ 
olution of ^ 22 meV. CeTes single crystals were grown 
by slow cooling of a binary melt [1^. Measurements 
were performed on freshly cleaved CeTes samples at 
T = 100 K in an ultrahigh vacuum chamber with a pres¬ 
sure of < 1 X 10“^^ torr. Phonon eigenmodes of the 
normal state of RTes were calculated from the dynami¬ 
cal matrix using the Quantum-ESPRESSO [33| package 
with ultrasoft pseudopotentials and the Perdew-Burke- 
Ernzerhof exchange-correlation functional, a kinetic 
energy cutoff of 500 eV, and 13x3x13 /c-points. Inci¬ 
dent pump fluences of 55 — 354 /iJ/cm^ correspond to 
absorbed energies of 16 — 103 meV per unit cell at the 
surface [35[. 2A decreases by no more than 15% and 
the system responds linearly [33 without melting of the 
charge order. This is in contrast to previous studies at 
higher excitation densities 0, 0, , where qualitative 

changes in the dispersion prevent decoupling the effects 
of coherent in-plane and out-of-plane motions. 

Fig 111(b) illustrates the cut through the ES where the 
trARPES data, which are shown in Eig. [TKc) and (d) for 
selected pump-probe delays, have been taken. Due to an 
asymmetric band dispersion and a slightly curved cut in 
the k = {kx^kz} plane the Eermi momentum kp does not 
exactly coincide with the maximum (minimum) position 
of the occupied (unoccupied) band. Eor simplicity, dis¬ 
persions along k are projected on kx- At a pump-probe 
delay of —330 fs the Te-5p band with the occupied part of 
the CDW gap is observed in equilibrium. At —25 fs, the 
leading edge of the pump pulse promotes electrons into 
unoccupied Te-5p states without yet changing the band 
dispersion significantly, which occurs at later delays. The 
magnitude of the full CDW gap 2A = 0.59(2) eV is di¬ 
rectly given by the difference of the maximum of the 
lower band at E — Ep = —0.16(1) eV and the minimum 
of the upper band at E — Ep = 0.43(1) eV. This is in 
agreement with previous results from optical reflectivity 
I^,[l^ while being smaller than extracted from ARPES 
21 and larger than estimated by tunneling experiments 

i3- 

To obtain a quantitative description of the electronic 
bands near Ep we consider a TB model of the in-plane 
Te-hpx/z orbitals [0,[2^ as sketched in Eig.[T]( a). Nesting 
of the metallic bare bands e{k) via q creates two branches 
which are separated by 2A 0: 


^1,2 (fc) = ± 


^k — q 


+ A2. (1) 
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FIG. 2. (a) Measured dispersion of the occupied Te-5p band 

for selected delays {markers) and transient TB fits {lines). 
Sketched effect of a transient rigid energy shift F^o(t) (b) and 
change of the CDW gap A(t) (c) on the band dispersion, (d) 
2A(t) and F^o(t) with smooth backgrounds {blaek lines) of a 
10^^ order polynomial indicated. 


e(k ) is approximated by a non-interacting TB model 
|24| of the in-plane Te-5p orbitals, see Fig. [T]^a). 

e{k) = - 2tpa cos{{kx - k^)al2) (2) 

- cos{{kx + kz)al2) - ^ Eq . 


The TB coupling parameters tpa- = —1.9 eV along 
and tpj^ = 0.35 eV perpendicular to the Te-chains 
(Fig . [Ha)) have been established by ARPES experiments 


22|, [2^. Ep = — 2tpcr sin (tt/S) and = 0.685 (27r/a) 


are constant. The in-plane lattice spacing is given by 
a = 4.34 A [l5|. The band dispersion is extracted by 
fitting a Gauss function to the spectra as function of 
k, and yield the data points in Fig. Hl^e). Eqs. (HD 
and (|2j) are combined and simultaneously fitted to the 
extracted dispersion above and below Ep at —25 fs as 
shown in Eig.ITJe). The rigid energy shift Eq = 0 is kept 
fixed and the only free fit parameter is A. This yields 
2A = 0.59(1) eV, which is identical to the value derived 
from the difference of upper and lower band energies at 
-25 fs. 


After having established that the TB model is well 
suited to describe the gapped dispersion near Ep, it is 
applied to the occupied band for all delays, as shown 
in Eig. |2fa), with A{t) and Eo{t) now being free time- 
dependent fit parameters. The low excitation regime al¬ 
lows for holding the TB parameters tpa and tpj^ constant. 
A{t) leads to momentum-dependent modulations in the 
spectra, and is naturally sensitive to the AM B as il¬ 
lustrated in Eig.ISfc). In contrast, Eo{t) is a momentum- 
independent energy shift of the whole band structure, as 
sketched in Eig.ISfb). As we will show, these parameters 
capture the two types of modes observed. 


Results for 2A(t) and Eo{t) are shown in Eig. ISfd). 
The magnitude of Eo{t) amounts to a few meV and is 
much smaller than the changes in 2A(t), which are on 
the order of several 10 meV. Within the pump pulse du¬ 
ration, 2 A drops by 15 % yet the charge order is only 
perturbed and not destroyed [35|. After ^ 3 ps, the sys¬ 
tem has thermalized at reduced gap size, which agrees 
with the timescale of the suppression of the structural 
order parameter satellite observed in time-resolved elec¬ 
tron diffraction 36, 
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FIG. 3. (a) Residuals of the polynomial background fit 

shown in Fig.lJJJd). (b) and (c) FT of (a). The vertical lines 
and markers at = 2.19(1) THz, Q 2 = 2.69(1) THz and 
Q 3 = 2.98(1) THz indicate mode frequencies and amplitudes 
independently obtained from fitting damped cosine functions 


The pronounced beating patterns of 2A{t) and Eo{t) 
exhibit an oscillatory response with similar yet distinct 
frequencies. For a more detailed analysis smooth back¬ 
grounds (10^^ order polynomials) are subtracted and the 
residuals 62A{t) and 6Eo{t) are plotted in Fig. |3fa). 
Analysis via Fourier transformation (ET) of 62A{t) re¬ 
veals two dominant frequencies at 1 = 2.2(1) THz and 
II2 = 2.7(1) THz, as shown in Eig. life). In contrast, the 
oscillatory response and particularly the beating pattern 
of 6Eo{t) differs from (52A(t). In addition to Hi and H2, 
the ET of 6Eo{t) exhibits a third pronounced peak at 
H3 = 3.0(1) THz, marked by a vertical arrow in Eig. 
|3](b). These findin gs a re robust regardless of the back¬ 
ground subtraction |35j. Fitting damped cosine functions 
to 62A{t) and 6Eo{t) [33 yields the same frequencies al¬ 
beit with higher resolution as indicated in Fig. |3l 

The mode at Hi is identified as a previously observed 
AM of RTes. As sketched in Fig.lU^a), coherent in-plane 
atomic motions modulate the lattice distortion Az and 
thus directly affect A. This assignment is supported by 
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comparing frequency to the soft mode observed in 
temperature-dependent transient optical reflectivity [?, 


my 

temperature-dependent Raman spectroscopy [l9|, l25| 


■26 1. and temperature- and fluence-dependent trARPES 

TOQ. 

The mode at Q 2 is missing at the zone center of the 
frozen phonon calculations, suggesting it originates either 
from a normal state optical or acoustic branch which is 
renormalized by coupling to the charge modulation at 
q. A flnite projection of the ionic motion on the lattice 
displacement Az associated with the CDW then natu¬ 
rally explains its presence in A{t) and the temperature- 
dependent softening in 26, [l^. We thus assign as a 
2nd AM. 


The mode at Qs is assigned to an out-of-plane Ai^ 
phonon, which is not related to the CDW by compari¬ 
son to frozen phonon calculations and transient reflectiv¬ 
ity [l^. The calculated displacement of the Te-atoms is 
sketched in Fig. SJ^c). The response of the Te-5p band is 
remarkably mode-dependent; Qs mainly affects Eq, and 
has a much smaller inffuence on A even though Qs in¬ 
volves the motion of the Te atoms, which form the CDW 
bands. However, the out-of-plane motion couples only 
weakly to the charge order, with correspondingly small 
inffuence on the order parameter A. 

To further investigate this explanation for mode- 
dependent coupling we plot FT amplitudes for all three 
modes in Fig. |4] as function of ffuence and momentum. 
These are derived from an FT analysis of the oscillatory 
binding energies, which have been extracted from tran¬ 
sient spectra as function of k [flU. The inffuence on the 
dispersion is illustrated in Fig. m^a) and (c). The AM 
at fti is strongly ffuence- and momentum-dependent as 
expected for a reduction of A. The AM at Q 2 exhibits 
a qualitatively similar behavior albeit of smaller magni¬ 
tude. This is consistent with a smaller coupling between 
the phonon with the charge density or a smaller projec¬ 
tion of the ionic motion on the lattice distortion A2;. 


In contrast, the out-of-plane Ai^ phonon at Qs causes 
a strikingly /^-independent shift, which is only weakly 
dependent on ffuence. This suggests a small ffuence de¬ 
pendent atomic displacement /i(t) and a constant defor¬ 
mation potential D^p(k) that explains the rigid energy 
shift via Fo(t) ^ D^p(k)ii(t) [H, [s^. We confirm the 
approximation of a k-independent deformation potential 
by DFT calculations for atomic displacements /i of up to 
0.13 A, shown in Ref. [sHj . 

In summary, trARPES measurements on CeTes reveal 
a CDW gap size of 2A = 0.59 eV. Coherent excitations 
of two AMs at Qi = 2.2 THz and at Q 2 = 2.7 THz and a 
generic Ai^ phonon at = 3 THz are identified by fits 
of a time-dependent model dispersion and comparison 
to normal state frozen phonon calculations. The out-of- 
plane Ai^ phonon leads to small rigid shifts of the entire 
Te-5p band but does not affect A. In contrast, and 
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FIG. 4. Mode-dependent modulations of the Te-5p bands 
(markers) by the three dominant modes is derived from 
momentum-dependent FT amplitudes for deposited energies 
of 28 — 103 meV per unit cell. Solid lines are obtained from 
the corresponding TB fits of the transient Te-5p bands 
(a) The AM at = 2.2 THz modulates the in-plane lat¬ 
tice displacement Az and leads to a strongly momentum- and 
fluence-dependent response, (b) The AM at D 2 = 2.7(1) THz 
exhibits a weaker momentum- and fluence-dependence, which 
indicates a weaker influence on the in-plane CDW. (c) The 
out-of-plane Ai^ mode at D 3 = 3 THz only shifts the band 
rigidly as it does not couple to the in-plane charge order. 


O2 result in pronounced changes of A. We point out that 
our /^-dependent analysis combined with low excitation 
densities is key to identifying and assigning AMs. We 
have generalized our approach and introduced a model- 
free evaluation (Fig|4)), which can be applied to materials 
where a detailed description of the band dispersion may 
not be available. More specifically, our approach can 
be transfered to charge-ordering in striped cuprates 0, 
which has identical symmetry to RTes [M |. 
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